Purcell factor for point— like dipolar emitter coupling to 2D-plasmonic waveguides 
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We theoretically investigate the spontaneous emission of a point-like dipolar emitter located near 
a two-dimensional (2D) plasmonic waveguide of arbitrary form. We invoke an explicite link with the 
density of modes of the waveguide describing the electromagnetic channels into which the emitter 
can couple. We obtain a closed form expression for the coupling to propagative plasmon, extending 
thus the Purcell factor to plasmonic configurations. Radiative and non-radiative contributions to 
the spontaneous emission are also discussed in details. 

PACS numbers: 42.50.Pq, 42.50.Nn, 73.20.Mf , 42.82.-m, 32.50.+d 



In 1946, Purcell demonstrated that spontaneous emis- 
sion of a quantum emitter is modified when located in- 
side a cavity [1]. A critical parameter is the ratio Q/K//, 
where Q and Vef / refer to the cavity mode quality factor 
and effective volume, respectively. In the weak coupling 
regime, the Purcell factor Fp, quantifies the emission rate 
7 inside the cavity compared its free-space value 70 
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where A is the emission wavelength and ni the cavity 
optical index. When Q/Veff is high enough, strong cou- 
pling regime occurs with reversible energy exchange be- 
tween the emitter and the cavity mode (Rabi oscillations) 
[2]. The design of cavities maximizing this ratio in order 
to control spontaneous emission is extremelly challeng- 
ing. There is however a trade-off between Q factor and 
effective volume. On one side, ultra high Q (~ 10^) are 
obtained in microcavities but with large effective volume 
(^ 10^ jum^). On the other side, diffraction limited mode 
volume [Veff ~ (A/ni)^] are achieved in photonic crys- 
tals but at the price of weaker quality factors {Q ~ 10^). 
Moreover, it is sometimes preferable to optimize Q I Veff 
but keeping a reasonable Q factor in order to efficiently 
extract the signal from the cavity. Additionally, the emit- 
ter spectrum can be large at ambiant temperature and 
better coupling is expected with low Q cavities [3] (i.-e. 
matching cavity and emitter impedances Q). 

In this context, it has been proposed to replace the 
cavity (polariton) mode by a surface plasmon polariton 
(SPP) sustained by metallic structures as an alternative 
to cavity quantum electrodynamics [HI, Q . SPP can have 
extremelly reduced effective volume, insuring high cou- 
pling rate with quantum emitters, albeit a poor quality 
factor (Q ~ 100 [7]). Particularly, coupling an emit- 
ter to a plasmonic wire shed new light on manipulating 
single photon source at a strongly subwavelength scale, 
with applications for quantum information processing [sj . 
Others promissing applications deal with the realization 
of integrated plasmonic amplifier (qUTi) . Highly resolved 
surface spectroscopy was also pointed out based either 



on the antenna effect [l2j or coupling dipolar emission to 
an optical fiber via a plasmonic structure [ij, ll^ . 




FIG. 1: Practice models. A dipolar emitter p is located at 
distance d of an infinite silver cylinder of circular (a) or pen- 
tagonal (b) cross-section, c) The dipolar emitter is located 
in a substrate- wire gap. 



In this work, we present an original approach for cal- 
culating rigorously the coupling of dipolar emitter to 2D 
plasmonic waveguides of arbitrary profile. We achieve 
a closed form expression for the coupling rate into the 
guided SPP. We also investigate the radiative and non 
radiative channels. In particular, the contribution of the 
plasmon, difficult to estimate otherwise [H, [l5|, is clearly 
established. Our method is general and treat equivalently 
bound and leaky waveguides of arbitrary cross-section, 
possibly on a substrate (Fig.[T]). 

According to Fermi's golden rule, coupling of a quan- 
tum emitter to a continuum of modes is governed by the 
(3D) local density of states (3D-LD0S) 
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where pu(r,cj) is the local density of modes, projected 
along the direction of the dipolar transition moment 
p = pu (partial LDOS) 0. r is the emitter location 
and 00 its emission frequency. To characterize the cou- 
pling independantly of the emitter properties, we intro- 
duce the normalized quantity 7(r)/7o = pu(r, cj)/p2i(^) 
where Pu(^) = cj^/Gtt^c^ is the free-space partial LDOS. 
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Since we are interested in 2D waveguide, the main idea 
is to work on the density of modes associated with the 
guide (bound and radiation modes). For this purpose, we 
now estabhsh a relationship between 2D and 3D LDOS 
by introducing Green's dyad formahsm. First, the 3D- 
LDOS is related to the 3D Green's tensor G of the system 
(Im and Tr refer to the imaginary part and trace) [l7| 
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In presence of an infinitely long (2D) structure, the 3D- 
Green's tensor is expressed by a Fourier transform 
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Then, we obtain the 3D-LD0S as a function of 2D- 
Green's dyad 
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Equation ^ obviously reproduces the 3D-LD0S in 
a homogeneous medium of index ni. Limiting 
the integration range to radiative waves, and since 
-:^/mTrG^f^(r||,r||,/c^) = cj/27rc^ in a homoge- 
neous medium, we obtain, as expected, po(r) = 
^/-nX^^^ cj/27rc2 = niUJ^/27T^c^, The quantity 

— ;^/mTrG^^(r|| , ry , /c^) is generally referred as 2D- 
LDOS by analogy with 3D-LD0S expression §1 [18]. It 
is a key quantity to understand spatially and spectrally 
resolved electron energy loss spectroscopy [19]. Equation 
([5]) makes then a direct link between 2D and 3D LDOS. 
We however consider a slightly different definition, more 
appropriate to describe a density of guided modes ^] 

2k 

The 2D Green's dyad is separated in two contributions 
q2D ^ q2D ^ ^q2D ^^^^^ q2D 2D-Green's 

dyad without the waveguide and AG^^ is the guide 
contribution. This formulation separates the reference 
system (multilayer substrate, homogeneous background, 
. . . ) from the guiding structure. It comes, with e^g/ the 
dielectric constant of the reference system. 



p2^(r|| , k,) = '-ImTr e{ru )G^''{ru , r,, , k,) . (6) 
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Ap'^ = -^/mTr e(r||)AG^^(r||,r||,fe,). 

This wording separates the continuum of modes of the 
reference system from the waveguide density of 

modes Ap^^ . The partial 2D-LD0S is finally 
2k 

Apl"" (r II , A:, ) = --^ImTr ) [u • AG^^ (r y , r y , fc,) • u] 

(8) 
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FIG. 2: (Color online) a) 2D radial LDOS variation as a 
function of kz at two distances to the nanowire of Fig. [T^). 
b) Logscale over the high momentum range. R — 20 nm, 
62 — —50 + 3.85i, A = 1 /jm and a — 2. 



Figure [2] represents the radial 2D-LD0S Apl^ik^) for 
the benchmark model defined in Fig. [1^). 2D-Green's 
dyad has been numerically evaluated by applying a mesh- 
ing on the waveguide cross-section [20]. The main con- 
tribution is the Lorentzian variation peaked at the effec- 
tive index of the guided SPP ng// = kspp/ko = 2.28, 
and with a full width at half maximum inversely propor- 
tional to the mode propagation length Lgpp = 1.2 pm (in- 
set). For kz < ni/co, the 2D-LD0S describes scattering 
events and contributes to radiative rate jrad- Finally, for 
kz > ni/co, LDOS takes part to the non-radiative decay 
rate jnr- Indeed, the plasmon is dissipated by thermal 
losses. Moreover, for very short distances, the 2D-LD0S 
spectrum extends over very large values of kz (Fig. [2)3). 
This behaviour is typical for non-radiative transfer by 
electron- hole pairs creation in the metal [21]. 

The coupling rate into the propagative SPP is obtained 
using equations (|3|5|8p and keeping only the plasmon 
contribution by limiting the integration of Eq. (|5]) to kz 
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FIG. 3: (Color online) Variation of the rates as a function of 
distance to the silver nanowire for a radial dipole. a) Cou- 
pling rate into SPP obtained using i) our approach based on 
2D-LD0S formulation, including losses ii) exact lossless case 
and iii) quasi- static approximation, b) Radiation rate calcu- 
lated using 2D-LD0S formulation (solid line) or quasi-static 
approximation (dotted line), c) Comparison of the plasmon 
rate jpi with the total non radiative rate jnr- 



corresponding to the SPP resonance. This is strongly 
simplified by the Lorentzian shape of the resonance and 
leads to the closed form expression 



rate into the plasmonic channel as a function of distance 
to the silver nanowire obtained using i) closed form ex- 
pression (j9j), ii) exact expression for a lossless plasmonic 
waveguide (p!Q]) and iii) a quasi-static approximation 

Quite surprinsingly, although the exact expression ne- 
glects dissipation, we obtain an excellent agreement with 
our expression that correctly accounts for losses. In for- 
mula (j9j) the ratio Ap'^jLgpp is proportional to the num- 
ber of guided modes [20| so that it does not depends on 
the losses. When losses tends towards zero, Lspp oo 
and Ap^^ ^ oo at resonance so that Ap^^ /Lgpp remains 
constant (Dirac distribution). Equivalently, this simply 
reveals that the emitter couples to the guided mode, no 
matter if the energy is dissipated by losses during prop- 
agation or propagates to infinity. 

We now turn on the radiative decay rate associated 
with the 2D-LD0S in the interval [—niko : niko]. 
We compare in Fig. [SJd) our numerical simulation with 
quasi-static approximation derived in Ref. [HI, for 
the nanowire. The quasi-static approximation underes- 
timates the radiative contribution to the coupling rate 
since it only considers the cylindrical dipole mode. 

Finally, the non-radiative decay rate jnr is deter- 
mined from 2D-LD0S calculated on the evanescent do- 
main \kz\ > niko which includes all the non radia- 
tive mechanisms: Joule losses during plasmon propa- 
gation and electron-hole pairs creation into the metal. 
Figure [3fc) represents the plasmon and total non ra- 
diative rates. The non-radiative rate diverges close to 
the wire surface whereas plasmon contribution remains 
finite. For large separation distances, the plasmon is 
the only contribution to the non radiative rate. We 
achieve an optimal coupling efficiency into the guided 
SPP, f3 = -fpi/{jrad + Inr) = 83%, at d = 20 nm. 
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This important result describes the emitter coupling rate 
to a 2D waveguide of arbitrary cross section. It is ex- 
pressed as the overlap between the dipolar emission and 
the guided mode profile {Ap^) divided by the mode 
propagation length in the longitudinal direction. This de- 
fines the 3D Purcell factor for a 2D geometry. Although 
presented for plasmonic waveguide, the demonstration 
remains valid for any 2D configuration (plasmonic cav- 
ity or waveguide pjj], metal coated or dielectric 
[22| nanofiber, ...). In order to validate this expression, 
we now compare it to the exact expression obtained by 
considering coupling to a lossless waveguide [H, : 



jp^ ^ 37rcE^{d)[E^{d)Y 
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FIG. 4: Coupling rate to guided SPP calculated near a cylin- 
drical wire of circular (solid line) or pentagonal (dotted line) 
cross-section (R = 20nm) . The modes profiles are shown. 



where (E, H) is the electromagnetic field associated with 
the guided SPP. In Fig. [3^), we compare the coupling 



So far, we considered a silver circular nanowire em- 
bedded in a homogeneous background to illustrate and 
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validate our method. In the fohowing, we investigate the 
two complex geometries depicted on Fig. [l]3,c). Figure 
m presents the coupling rate into the SPP supported by 
a penta-twinned crystalline nanowire recently character- 
ized [in. At short distances, the coupling rate into the 
guided SPP is strongly enhanced as compared to coupling 
to a circular wire of similar dimensions. This is due to 
the strong mode confinement near the wire corners as 
revealed by the mode profile. 
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FIG. 5: (Color online) Different contributions to the decay 
rates for a 100 nm diameter silver wire 50 nm above a glass 
substrate (es = 2.25). Superstrate is air (ei = 1). 

Experimental configurations generally concern struc- 
tures deposited on a substrate. For high index substrate, 
the otherwise bound mode becomes leaky. Note that 
usual expression (p!Q|) is then practically unenforceable 
due to difficulty of normalizing the mode. Differently, 
expression (j9j), derivated in this work, is easily used even 
in such a situation. Moreover, in case of leaky mode, 
it is even more difficult to properly distinguish radiative 
and non radiative contributions to the coupling rate, as 
compared to the bound mode situation treated above. In- 
deed, the guided plasmon contributes to both the radia- 
tive rate (leaky part) and non radiative transfer (intrinsic 
losses). This difficulty is easily overcome using the 2D- 
LDOS formalism. The propagation length can be written 
Lspp = {^racf + ^nr(ui)~'^ whcrc the radiative and non 
radiative rates have been introduced. As an example, we 
consider a 100 nm silver wire 50 nm above a glass sub- 
strate. We calculate an effective index ng// = 1.28, be- 
low the substrate optical index, indicating a leaky mode. 
Its propagation length is Lspp = 1.2 /im = 
with F*^^^ = 0.083/im~^. The leakage rate is evaluated 
by cancelling the metal losses (/m(e2) = 0). We ob- 
tain F;^^^ = 0.073/im~^. Figure [5] shows the interplay 
between the various contributions to the decay rate for 
an emitter placed in the wire-substrate gap. The radia- 
tive rate -frad = Iscatt + IpiMak IS the sum of the scat- 



tering and leakage channels and the non radiative rate 
iNR = Jpi,NR -^le-h Originates from plasmon losses and 
electron-hole pairs creation. Except for short distances, 
the main decay channel is the plasmon decoupling into 
the substrate. We obtain a maximum decoupling emis- 
sion into the substrate j3 = ^pi^ieak/l — '^0% for an emit- 
ter centered in the gap [d = 25 nm) [26^]. 

To conclude, we derive an explicit expression for the 
coupling rate between a point-like quantum emitter and 
a 2D plasmonic waveguide. We define the coupling Pur- 
cell factor into the plasmon channel whereas the radia- 
tive and non radiative rates are numerically investigated. 
This method clearly reveals the physics underlying the 
complex mechanisms of spontaneous emission coupled to 
a plasmonic guide (scattering, leakage, electron-hole pairs 
creation, SPP excitation). 
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